Designing bound states in a band as a model for a quantum network 
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We provide a model of a one dimensional quantum network, in the framework of a lattice using 
Von Neumann and Wigner's idea of bound states in a continuum. The localized states acting as 
qubits are created by a controlled deformation of a periodic potential. These wave functions lie at 
the band edges and are defects in a lattice. We propose that these defect states, with atoms trapped 
in them, can be realized in an optical lattice and can act as a model for a quantum network. 
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In classical computation, bits and 1 are used to represent data. Logic gates are employed for computation and 
data manipulation. At present, considerable work is being done to develop quantum computers. It is expected that 
, they will increase the speed and efficiency of computation [Q compared to a classical computer. The basic building 
blocks in quantum computing are known as quantum bits or qubits and are represented by the eigenstates of a system. 
For example, qubits are described by the up and down states of spin-half particles, represented by |0) and |1). The 
unitary operators play the role of logic gates and we can perform various operations similar to classical computing. 
The major areas of interest in the field of quantum computation are the development of models, having states which 
can be used as qubits and the construction of unitary operators to manipulate these qubits. It is important that these 
states should be well isolated from the outside environment to avoid environmental decohcrence. 

In this paper, we propose a theoretical model for preparing qubits i.e., |0 > and |1 >, using periodic potentials. 
• These potentials are characterized by a band spectrum, which contains bands of allowed energies interspersed with 
forbidden energy gaps. The wave functions are not integrable and they extend to spatial infinity. Hence, they cannot 
r-j 1 be envisaged as qubits. We will explicitly demonstrate that a periodic potential can be appropriately deformed to 
fS |' accommodate bound states in its band spectrum. This deformation can be treated as the perturbation of the original 
potential. The effect of which is the creation of the localized states which can be perceived as defects in the lattice. 
Atoms trapped in these defects can be used to describe qubits. This treatment of the trapped atoms in optical 
lattices, as qubits in a quantum network, is similar to the quantum computation model proposed by Angelakis et. 
ai.,P]. Here, they have used photonic crystals to confine photons in the defect states created inside the band gap, 
to represent qubits. In quantum networking atoms are trapped in the nodes and are used to store information. We 
propose that an array of such localized states in the deformed optical lattice, with atoms trapped in them, can be 
k> j used as a quantum network [00]. 

At present, various optical lattices are routinely realized in the laboratories and atoms can be trapped in a potential 
well with relative ease. For our model, we consider one of the well - known family of periodic potentials namely, the 
Lame potentials E] , 

V(x) — j(j + l)msn 2 (x,m). (1) 

The function sn (x,m) is the Jacobi elliptic function [7] [SI with elliptic modulus < m < 1. These potentials are 
exactly solvable for integer values of j and for a given j there are (2j + 1) band-edges. Both the form of these 
band-edge solutions and explicit solutions for smaller values of j are given in El El QDI • The Lame potential has been 
proposed as a model for quasi 1 - d confinement of Bosc - Einstein condensates (BEC) in a standing light wave \1'2\ . 
Possible application of BECs for quantum computation are currently being explored and the bound states created by 
deforming a potential like the Lame potential may have useful implications. 

The method used to construct these bound states is the same as that used by Pappademos et. at, to construct 
quantum mechanical bound states in a classical continuous energy spectrum. Such states were first discovered by 
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Von Neumann and Wigner ^3] an d many such examples were found later ^]-^]. It became clear that such states 
arose due to the delicate interplay between the rate at which the oscillatory potential falls off and the time taken by 
the various maxima of the potential to interact and create a bound state. Existence of such states was reported by 
Capasso et. ai, in semiconductor hetrostructures |17j . 

SUSYQM ^| EH was applied to generate bound states in the continuum of the energy spectrum for a spherically 
symmetric potential |20l I21j . We apply the same technique to periodic potentials and construct bound states in their 
band spectrum, with the same energies as the band-edge energies. The solutions thus obtained are located at the 
band-edges and not inside the forbidden energy gap as discussed in [21 03 HID ■ 

In the next section, we give an overview of SUSYQM and briefly describe the steps involved in constructing a 
bound state in the continuum (BIC). In section III, we use the same procedure to construct bound states for the 
Lame potential. In the last section, we present our conclusions. 



II. SUPERSYMMETRIC QUANTUM MECHANICS 



For a given 1 - d potential V(x) with eigenfunctions u n (x) and eigenvalues E n (n = 0, 1,2...), we can generate a new 
family of potentials V(x;X), isospectral to V{x). The parameter A is used to label the potentials in the isospectral 
family and takes values lying in the range A > and A < — 1. Setting h — 2m = 1 we can write the superpotential as, 



W{x) = -iS. (2) 

Uq 

The original potential V(x) can be expressed in terms of W(x) as, 

V{x) = W 2 {x) - W'{x). (3) 

Its isospectral partner V(x; A) is given by 

V{x;X) = W 2 {x) + W\x). (4) 

Let u n (x) be the eigenfunctions of V(x; A) (n — 1, 2...). We have u n {x) = Au n (x) where A = d/dx + W(x). (Hence, 
A' = —d/dx + W(xj). Note that uq(x), the ground state cannot be obtained in this manner since Auo(x) — 0. Thus, 
V(x; A) is isospectral to V{x), except that its spectrum does not contain uq(x). Hence, to introduce the ground state 
into its spectrum and form a complete set of eigenstates, we need to find the most general superpotential W(x), so 
that 

V(x;X) = W 2 {x) + W'(x). (5) 

We can now show that 

W(x) = W(x) + 4- Woix) + A) (6) 
dx 

where 

7 = / ul{y)dy. (7) 
Jo 

Thus, in this process of reinstating the ground state, we obtain the expression for the potential V(x; A) in terms of Jo 
as, 

V{x- A) = V(x) 2[ln(J + A)]" = V{x) - ^ + . (8) 

I + A (J + A)^ 

Note that V(x; A) is V(x) plus a perturbative term and A can be take as the perturbative parameter. The perturbed 
potential is isospectral to the old potential and its eigenstates are given below in Eqs. © and (|10|) . 

In standard SUSYQM, isospectral families of potentials which allow only bound states have been constructed, 
and uq(x) was taken to be the ground state. In |21|. the above method was generalized to the case where the 
potentials have a continuous energy spectrum and uq(x) has been taken to be any non-singular eigenstate of V(x). 
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This procedure was used to construct BIC for the spherically symmetric potential and the results were stated as a 
theorem. We refer the reader to [2D] and the references therein for further details and merely reproduce the result here. 

Let uq{x) and u\{x) be any two nonsingular solutions of the Schrddinger equation for the potential V(x) correspond- 
ing to arbitrarily selected energies Eq and E\ respectively. Construct a new potential V(x; A) as prescribed by Eq. Jig) . 
Then, the two functions 

JO + A 

and 

ui(x) = (Ei -E )ui(x) + u (x)W r (u (x),ui(x)) (10) 

are solutions of the Schrddinger equation for the new potential V(x;X), corresponding to the same energies Eq and 
Ei. Here W t {uq(x), Ui(x)) is the Wronskian. 

Note that the original potential V(x) had no integrable solutions but the new potential has one square-integrable 
solution Cto (x) , with the rest being non - integrable. The creation of the bound state can be elucidated by the fact 
that Io in Eq. J7J diverges owing to the non-integrability of uq. Hence, as Iq — > oo, uo(x) — ■> 0, resulting in a 
square-integrable wave function in the continuum. 

We can create another bound state by using the non-normalizablc state iii(x) in place of uq(x) and deforming 

V(x; A), using the same procedure described above. We then obtain a potential V(x; A, Ai), isospectral to V(x) which 
has two bound states in the continuous spectrum with energies Eq and E\. The parameter Ai is a real number lying 
in the range Ai > and Ai < —1. 

The expressions for the new potential V(x) and the two square-integrable states uq(x) and ui(x), with energies Eq 
and Ei respectively are given by 

V{x) = V{x) 2[ln(7 x + Ai)]" = V{x) - i*j^L + 2 f (11) 

7i + Ai (Ji + Ai) 2 

uq(x) = (E - Ei)uo(x) + ui(x)W r (ui(x), u (x)) (12) 

and 
where 

h= f u 2 i(y)dy (14) 
Jo 

and W r (ui(x), uq(x)) is the Wronskian. In the next section, we apply the above technique to periodic potentials on 
the half-line and construct bound states in the band spectrum. For this purpose, we use the band-edge wave functions 
to deform the original periodic potential. The bound states thus created have the same band-edge energies. 

III. THE LAME POTENTIAL 

The Lame potential (Eq. (Q) with j = 2 is, 

V(x) = 6msn 2 (x, m). (15) 

It has two bands and a continuum. The expressions for the five band - edge wave functions and energies are given 
below [H| [l?fl 110) |27) with ipo(x) and ipi( x ) representing the lower and upper band - edge wave functions of the first 
band and so on. We have 

i/j (x) = 3m + 3- 6 -3msn 2 (x,m) , E = 2S-2m-2, (16) 
tpi(x) — cn (x, to) dn (x, m) , Ei=m + 1, (17) 
^{x) = dn (x, to) sn (x, to) , E 2 = 4771 + 1, (18) 
ip3 (x) = cn (x, to) sn (x, to) , E% — m + 4, (19) 
ipi(x) = 3TO + 3-35-3TOsn 2 (x,TO), £4 = 25 + 2to + 2, (20) 
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where, the functions cn(x,m) and dn(x, m) are the Jacobi elliptic functions with modulus parameter m. For con- 
structing the bound states, we examine only the half-line problem and hence consider the band-edge wave functions 
which vanish at the origin. In the above given expressions, only Eqs. IjlSJI and (|19fl which represent the lower and 
upper band edge wave functions of the second band, satisfy this condition, since sn (0, m) = 0. We follow the steps 
described in the previous section and construct two bound states with energies Ei and E%. The entire procedure is 
done numerically and we give the plots of the deformed potential and bound states thus obtained, in the sequel. 

We deform the potential given in Eq. i|15fl . using the expression for ^(x) in Eq. (|18|) . to obtain a one-parameter 
bound state solution which depends on the parameter A . For this purpose, we first plot Iq versus x in Fig. ^ 
As expected Iq turns out to be a diverging integral. Using Iq and Eqs. © and (|10|l . we plot the deformed 

potential V{x) and the deformed wave functions ip2(x) and ^{x). These are give in Figs. El El and E] respectively. 
For comparison, the original potential and wave functions are plotted in dotted line. 

It is clear from Figs. El and E] that ip2(x) is a normalizable state and ip3(x) is not normalizable. Thus, with this 
deformation we have obtained only one bound state. In order to construct two bound states, we deform V(x) with 
ip3(x) using Eqs l(TT|) - lfH|) . Plots of 1%, V{x), ^{x) and 4>z{x) versus x are given in Figs. El El and El respectively. 

From Figs[7|andEl it is clear that both ip2(x) and ipsix) are integrable. These states have energies E 2 and E 3 and 

the potential V(x) is isospectral to the original potential V(x). The deformed potential and the bound states depend 
on the parameters A and Ai and as we increases their values, the deformed states and the potential tend towards the 
corresponding original states and the potential. In Figs. El and we give the plots of the deformed wave functions 

■02 (x) and ips (x), for two different values of A and Ai- (Without loss of generality we have set A = Ai = 1 and 
A = Ai = 10 in these figures). 

IV. CONCLUSIONS 

We have shown that we can construct bound states in the band spectrum of a periodic potential using SUSYQM. 
It is clear from the above procedure that we can create a class of bound states by successively deforming the potential 
V(x), provided there exist band-edge wave functions, of the original potential, which satisfy the boundary condition 

iMo) = o. 

We can use these square-integrable states as qubits in quantum computation. The localized states with atoms 
trapped in them can be treated as the qubits |0) and |1). If it is possible to get n such defect states in the optical 
lattice we have an array of qubits which can be used as an optical network py. The perturbative parameters can 
be used to adjust the overlap of these trapped atoms and also to control the deformation of the periodic potential. 
Moreover, these localized states are in the band and hence, protected from the external influences. Since optical 
lattices are easy to create and manipulate in the laboratory, this can be a useful model for quantum networking |24l 
E3EH1- 

In conclusion, we have shown that it is possible to deform a periodic potential to accommodate localized states at 
the band edges. It is proposed that by trapping atoms in these states we can construct qubits and an array of n such 
qubits can be used as a quantum network. 
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FIG. 1: Plot of the diverging integral Io versus x. 
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FIG. 2: Deformed potential V(x) versus x, with A = 1. The dotted line represents the original potential V(x). 
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FIG. 3: Deformed wave function ip2{x) versus x, with A = 1. The dotted line represents the original band-edge wave function 

1p2 (x). 
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FIG. 4: Deformed wave function 1^3(0;) versus x, with A = 1. The dotted line represents the original band-edge wave function 

1p3 (x). 
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FIG. 5: Plot of the diverging integral 7i versus x. 
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FIG. 6: Deformed potential V(x) versus x, with Ai = 1. The dotted line represents the original potential V(x). 
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FIG. 7: Deformed wave function ip2(x) versus x, with Ai = 1. The dotted line represents the original band-edge wave function 

4>2(x). 
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FIG. 8: Deformed wave function %p3{x) versus x, with Ai = 1. The dotted line represents the original band-edge wave function 
ip 3 {x). 
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FIG. 9: Deformed wave function ip2{x) for Ai = A = 1 (dotted line) and for Ai = A = 10 (thick line). 
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FIG. 10: Deformed wave function ip3(x) for Ai = A = 1 (dotted line) and for Ai = A = 10 (thick line). 



